Monogenic functions in 5-dimensional spacetime used as first principle: gravitational 
dynamics, electromagnetism and quantum mechanics 
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Monogenic functions are functions of null vector derivative and are here analysed in the geometric 
algebra of 5-dimensional spacetime, G4,i, in order to derive several laws of fundamental physics. The 
paper introduces the working algebra and the definition of monogenic functions, showing that these 
generate two 4-dimensional spaces, one with Euclidean signature and the other one with Minkowski 
signature. The equivalence conditions between the two spaces are studied and relativistic dynamics, 
not entirely coincident with Einstein's general theory of relativity, is demonstrated. The monogenic 
condition is then shown to produce Maxwell's equations and electrodynamics both classical and 
quantized. 
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I. INTRODUCTION 

Our goal is to show how the important equations of 
physics, such as relativity equations and equations of 
quantum mechanics, can be put under the umbrella of 
a common mathematical approach^Q- We use geomet- 
ric algebra as the framework but introduce monogenic 
functions with their null derivatives in order to advance 
the concept. Furthermore, we clarify some previous work 
in this direction and identify the steps to take in order 
to complete this ambitious project. 

Since A. Einstein formulated dynamics in 4- 
dimensional spacetime, this space is recognized by the 
vast majority of physicists as being the best for formu- 
lating the laws of physics. However, mathematical con- 
siderations lead to several alternative 4-D spaces. For ex- 
ample, the 4-dimensional space called 4-D optics (4D0) 
is equivalent to the 4-D spacetime of the general theory 
of relativity (GTR) when the metric is static, and there- 
fore the geodesies of one space can be mapped one-to-one 
with those of the other. Then one can choose to work in 
the space that is more suitable. 

In the case of a central mass, we can examine how the 
Schwarzschild metric in GTR can be transposed to 4D0. 
The usual form of the metric is 



At' 



1 ]dt^ - [1 

X J \ X 



dx' 



(1) 



where m is the spherical mass and x is the radial coor- 
dinate, not the distance to the centre of the mass. This 
form is non-isotropic but a change of coordinates can be 
made that returns the expression to isotropic form (see 
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D'Inverno '3', section 14.7]): 

r= (x~m + \Jy^ ~ 2mx) /2; 
and the new form of the metric is 
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From this equation we immediately define two coeffi- 
cients, which are called refractive index coefficients, 
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We devote the first part of this paper to deriving 
them from a geometric algebra approach in a special 
5D space with null geodesies, thereby establishing that 
there is a 4D0 Euclidean metric space equivalent to 
the Schwarzschild metric space. We build upon previ- 
ous work by ourselves and by other authors about null 
geodesies, regarding the condition that all material par- 
ticles must follow null geodesies of 5D space: 

The implication of this for particles is 
clear: they should travel on null 5D geodesies. 
This idea has recently been taken up in the 
literature, and has a considerable future. It 
means that what we perceive as massive par- 
ticles in 4D are akin to photons in 5D.||J| 

Accordingly, particles moving on null 
paths in 5D [dS'^ — 0) will appear as mas- 
sive particles moving on timelike paths in 4D 
(ds^ > 0) ... 

We actually improve on these null displacement ideas 
by introducing the more fundamental monogenic condi- 
tion, deriving the former from the latter and establishing 
a common first principle. 
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II. SOME GEOMETRIC ALGEBRA 

Geometric algebra is not usually taught in university 
courses and its presence in the literature is scarce; good 
reference works are 0, S 0- We will concentrate on 
the algebra of 5-dimensional spacetime because this will 
be our main working space; this algebra incorporates as 
subalgebras those of the usual 3-dimensional Euclidean 
space, Euclidean 4-space and Minkowski spacetime. We 
begin with the simpler 5D flat space and progress to a 
5D spacetime of general curvature (see Appendix |0 for 
more details.) 

The geometric algebra G4^i of the hyperbolic 5- 
dimensional space we consider is generated by the co- 
ordinate frame of orthonormal basis vectors such that 

(^o)' = -l, 

{<^.? - 1, (5) 
-cT/s = 0, a ^ (3. 

Note that the English characters i, j, k range from 1 to 
4 while the Greek characters a, (3, 7 range from to 4. 
See the Appendix^for the complete notation convention 
used. 

Any two basis vectors can be multiplied, producing 
the new entity called a bivector. This bivector is the ge- 
ometric product or, quite simply, the product, and it is 
distributive. Similarly to the product of two basis vec- 
tors, the product of three different basis vectors produces 
a trivector and so forth up to the fivevector, because five 
is the dimension of space. 

We will simplify the notation for basis vector products 
using multiple indices, i.e. aacris = aafs- The algebra is 
32-dimensional and is spanned by the basis 

• 1 scalar, 1, 

• 5 vectors, ctq,, 

• 10 bivcctors (area), ctq^, 

• 10 trivectors (volume), ctq^^, 

• 5 tetravectors (4- volume), iaa, 

• 1 pseudoscalar (5- volume), i = (701234 ■ 
Several elements of this basis square to unity: 

(cT^f - i'yo^f = - (i^o)' = 1. (6) 

The remaining basis elements square to —1: 

(^0)' = (a„,)2 = = iia,r = i^ = -1. (7) 

Note that the pseudoscalar i commutes with all the other 
basis elements while being a square root of — 1 ; this makes 
it a very special element which can play the role of the 
scalar imaginary in complex algebra. 

In 5-dimensional spacetime of general curvature, 
spanned by 5 coordinate frame vectors (7ct, the indices 



follow the conventions set forth in Appendix^ We will 
also assume this spacetime to be a metric space whose 
metric tensor is given by 

gap^da-gp; (8) 

the double index is used with g to denote the inner prod- 
uct of frame vectors and not their geometric product. 
The space signature is ( — h + + -f), which amounts to 
saying that goo < and gu > 0. A reciprocal frame is 
defined by the condition 

g'^-gp^sy. (9) 

Defining g"^^ as the inverse of gap, the matrix product 
of the two must be the identity matrix; using Einstein's 
summation convention this is 

g'^-'gp^ = <5"0. (10) 

Using the definition ((SJ we have 

ig"''g^)-gf3 = S"p; (ii) 

comparmg with Eq. ® we determine ^"with 

g" = g^'^g-,- (12) 

If the coordinate frame vectors can be expressed as a 
linear combination of the orthonormed ones, we have 

ga = n^acrp, (13) 

where n^a is called the refractive index tensor or simply 
the refractive index; its 25 elements can vary from point 
to point as a function of the coordinates. Q When the 
refractive index is the identity, we have ga = era for the 
main or direct frame and g'^ = — ctq, = Ci for the 
reciprocal frame, so that Eq. ^ is verified. In this work 
we will not consider spaces of general curvature but only 
those satisfying condition l|13|l . 

The first use we will make of the reciprocal frame is for 
the definition of two derivative operators. In flat space 
we deflne the vector derivative 

V = <T"da. (14) 

It will be convenient, sometimes, to use vector derivatives 
in subspaces of 5D space; these will be denoted by an 
upper index before the V and the particular index used 
determines the subspace to which the derivative applies; 
For instance — a"^dm = cr^di + (7^82 + 0-^83. In 
5-dimensional space it will be useful to split the vector 
derivative into its time and 4-dimensional parts 

V = -aadt + a'd, = -a^dt + V. (15) 

The second derivative operator is the covariant deriva- 
tive, sometimes called the Dirac operator, and it is de- 
fined in the reciprocal frame g"' 

D = 5"5„. (16) 
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Taking into account the definition of the reciprocal frame 
©, we see that the covariant derivative is also a vector. 
In cases such as those we consider in this work, where 
there is a refractive index, it will be possible to define 
both derivatives in the same space. 

We define also second order differential operators, des- 
ignated Laplacian and covariant Laplacian respectively, 
resulting from the inner product of one derivative opera- 
tor by itself. The square of a vector is always a scalar and 
the vector derivative is no exception, so the Laplacian is 
a scalar operator, which consequently acts separately in 
each component of a multivector. For 4 -f 1 flat space it 
is 

= (17) 

One sees immediately that a 4-dimensional wave equation 
is obtained by zeroing the Laplacian of some function 

= (18) 

This procedure was used in Ref. for the derivation 
of special relativity and will be extended here to general 
curved spaces. 



III. THE MONOGENIC CONDITION 

There is a class of functions of great importance, called 
monogenic functions &\ , characterized by having null vec- 
tor derivative; a function ■0 is monogenic in flat space if 
and only if 

VV^ = 0. (19) 

A monogenic function is not usually a scalar and has by 
necessity null Laplacian, as can be seen by dotting Eq. 

with V on the left. We are then led to Eq. ((TH)l . 
which can also be written as 

du^ = doo^. (20) 

i 

This relation can be recognized as a wave equation in the 
4-dimensional space spanned by the ai which will accept 
plane wave type solutions of the general form 

V> = 0oe'(P""°+'\ (21) 

where TpQ is an amplitude whose characteristics we shall 
not discuss for now, (5 is a phase angle and Pa are con- 
stants such that 

J2ip^f~ipor^0. (22) 
i 

When working in curved spaces the monogenic condi- 
tion is naturally modified, replacing the vector derivative 



V with the covariant derivative D. A generalized mono- 
genic function is then a function that verifies the equation 

B^j = 0. (23) 

Similarly to what happens in fiat space, the covariant 
Laplacian is a scalar and a monogenic function must ver- 
ify the second order differential equation 

= 0. (24) 

It is possible to write a general expression for the covari- 
ant Laplacian in terms of the metric tensor components 
(see 1^ Section 2.11]) but we will consider only situations 
where that complete general expression is not needed. 

When Eq. is multiplied on the left by D, we are 
applying second derivatives to the function, but we are 
simultaneously applying first order derivatives to the re- 
ciprocal frame vectors present in the definition of D itself. 
We can simplify the calculations if the variations of the 
frame vectors are taken to be much slower than those of 
function ip so that frame vector derivatives can be ne- 
glected. With this approximation, the covariant Lapla- 
cian becomes = g°''^daf3 and Eq. ^ can be written 

g^^do^li^j = 0. (25) 

This equation can have a solution of the type given by 
Eq. (|^ if again the derivatives ofpa are neglected. This 
approximation is usually of the same order as the former 
one and should not be seen as a second restriction. In- 
serting Eq. (|21|l one sees that it is a solution if 

= 0. (26) 

This equation means that the square of vector p — g°'pa 
is zero, that is, p is a vector of zero length and is called a 
null vector or nilpotent. Vector p is the momentum vector 
and should not be confused with 4-dimensional conjugate 
momentum vectors defined below. 



IV. EQUIVALENCE BETWEEN 4DO AND GTR 
SPACES 

By setting the argument of constant in Eq. (|21|l and 
differentiating we can get the differential equation 

Pada;" = 0. (27) 

The Ihs can equivalently be written as the inner product 
of the two vectors p • dx = 0, where dx = gpdx^ is a 
general 5D elementary displacement. In 5D hyperbolic 
space the inner product of two vectors can be null when 
the vectors are perpendicular but also when the two vec- 
tors are null. Since we have established that p is a null 
vector, Eq. (|27() can be satisfied either by da; normal to 
p or by (dx)^ =0. In the former case the condition de- 
scribes a 3-volume called wavefront and in the latter case 
it describes the wave motion. Notice that the wavefronts 
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are not surfaces but volumes, because we are working 
with 4-dimensional waves. 

The condition describing 4D wave motion can be ex- 
panded as 



dap 



dx"dx^ = 0. 



(28) 



This condition effectively reduces the spatial dimension 
to four but the resulting space is non-metric because all 
displacements have zero length. We will remove this dif- 
ficulty by considering two special cases. First let us as- 
sume that vector is normal to the other frame vectors 
so that all goi factors are zeroed; condition l|28|) becomes 



(29) 



All the terms in this equation are scalars and we are 
allowed to rewrite it with (da;°)^ in the Ihs 



(da;0)2 =-^dxMx^ 
500 



(30) 



We could have arrived at the same result by defining a 
4-dimensional displacement vector 



da:% 



-1 



gidx' 



(31) 



and then squaring it to evaluate its length; u is a unit 
vector called velocity because its definition is similar to 
the usual definition of 3-dimensional velocity; its compo- 
nents are 



Vi = 



dx' 
d^' 



(32) 



Being unitary, the velocity can be obtained by a rotation 
of the cr4 frame vector 



V — RaAR. 



(33) 



The rotation angle is a measure of the 3-dimensional ve- 
locity component. A null angle corresponds to v directed 
along tT4 and null 3D component, while a 7r/2 angle cor- 
responds to the maximum possible 3D component. The 
idea that physical velocity can be seen as the 3D compo- 
nent of a unitary 4D vector has been explored in several 
papers but see 10]. 

Equation H31|) projects the original 5-dimensional space 
into a space with 4 dimensions, with Euclidean signature, 
where an elementary displacement is given by the varia- 
tion of coordinate x^ . In the particular case where g^ — 
(To the displacement vector simplifies to dx^v = gidx"^ 
and we can see clearly that the signature is Euclidean 
because the four gi have positive norm. Although it has 
not been mentioned, we have assumed that none of the 
frame vectors is a function of coordinate x". 

Returning to Eq. (|28|l we can now impose the condition 
that 54 is normal to the other frame vectors in order to 
isolate (dx"')^ instead of (dx")^, as we did before; 



(dx 



4\2 



544 



dx^dx" 



(34) 



We have now projected onto 4-dimensional space with 

signature (-1 ), known as Minkowski signature. In 

order to check this consider again the special case with 
go — (To and the equation becomes 



1 



9r, 



(dx4)2 = — (dx0)2 - 2IIiIidx"dx"; 

544 544 



(35) 



the diagonal elements gu are necessarily positive, which 
allows a verification of Minkowski signature. Contrary 
to what happened in the previous case, we cannot now 
obtain (dx^)^ by squaring a vector but we can do it by 
consideration of the bivector 



dx V ^ 



1 



5445 



44 



5p5 



'^dx''. 



(36) 



All the products gfj,g^ are bivectors because we imposed 
54 to be normal to the other frame vectors. When (dx'')^ 
is evaluated by an inner product we notice that go5'^ 
has positive square while the three gm9^ have negative 
square, ensuring that a Minkowski signature is obtained. 
Naturally we have to impose the condition that none of 
the frame vectors depends on x^ . Bivector i/ is such that 
v — vv = 1 and it can be obtained by a Lorentz trans- 
formation of bivector cro4. 



v = TCT04T. 



(37) 



where T is of the form T — exp(_B) and i? is a bivector 
whose plane is normal to a 4. Note that T is a pure ro- 
tation when the bivector plane is normal to both (Tq and 
0-4. 

In special relativity it is usual to work in a space 
spanned by an orthonormed frame of vectors 7^ such 
that (70)^ = 1 and (7m)^ = —1, producing the desired 
Minkowski signature The geometric algebra of this 
space is isomorphic to the even sub-algebra of G4_i and 
so the area element dx^i^ H36I) can be reformulated as a 
vector called relativistic 4-velocity. 

Equations H30|l and (|34|l define two alternative 4- 
dimensional spaces, those oi 4-dimensional optics (4-DO), 
with metric tensor —gij/goo and general theory of rela- 
tivity (GTR) with metric tensor —g^u/ga, respectively; 
in the former x*' is an afhne parameter while in the lat- 
ter it is x'* that takes such role. In fact Eq. H34II only 
covers the spacelike part of GTR space, because (dx'^)^ 
is necessarily non-negative. Naturally there is the limi- 
tation that the frame vectors are independent of both x'' 
and x"*, equivalent to imposing a static metric, and also 
that goi = 5^4 = 0. Provided the metric is static, the 
geodesies of 4D0 can be mapped one-to-one with space- 
like geodesies of GTR and we can choose to work on the 
space that best suits us for free fall dynamics. For a phys- 
ical interpretation of geometric relations it will frequently 
be convenient to assign new designations to the 5D co- 
ordinates that acquire the role of affine parameter in the 
null subspace. We will then make the assignments x^ = t 
and x^ = T. Total derivatives with respect to these co- 
ordinates will also receive a special notation: df /dt = f 
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and d//dr = /. Special units conventions used in this 
paper are detailed in appendix IBl 

Unless otherwise specified, we will assume that the 
frame vector associated with coordinate cc" is unitary and 
normal to all the others, that is go = ctq and goi = 0. 
Recalling from Eq. (|30|l . these conditions allow the def- 
inition of 4D0 space with metric tensor gij. Although 
we could try a more general approach, we would loose 
the possibility of interpreting time as a line element and 
this, as we shall see, provides very interesting and novel 
interpretations of physics equations. In many cases it 
is also true that is normal to the other frame vec- 
tors and we have seen that in those cases we can make 
metric conversions between GTR and 4D0; it will be in- 
teresting, however, to examine one or two situations with 
non-normal 94 and so we leave this possibility open. 

For the moment we will concentrate on isotropic space, 
characterized by orthogonal refractive index vectors gi 
whose norm can change with coordinates but is the same 
for all vectors. Normally we relax this condition by ac- 
cepting that the three g^ must have equal norm but g^ 
can be different. The reason for this relaxed isotropy is 
found in the parallel we make with physics by assign- 
ing dimensions 1 to 3 to physical space. Isotropy in a 
physical sense need only be concerned with these dimen- 
sions and ignores what happens with dimension 4. We 
will therefore characterize an isotropic space by the re- 
fractive index frame go — ctq, gm ~ n-rcrm, 54 ~ n^a^. 
Indeed we could also accept a non-orthogonal 94 within 
the relaxed isotropy concept but we will not do so for the 
moment. 

Equation H30|l can now be written in terms of the 
isotropic refractive indices as 



df = K)2 ^(dx'")2 + (n^dTf 



(38) 



Spherically symmetric static metrics play a special role; 
this means that the refractive index can be expressed as 
functions of r if we adopt spherical coordinates. The 
previous equation then becomes 



dt^ ^ {urf [dr^ + r^{dQ^ + sin^ M^?^)] 
+ {n^dTf. 



to Schwarzschild metric, allowing 4D0 to be used as an 
alternative to GTR. Recalling that we derived trajecto- 
ries from solutions H21|l of a 4-dimensional wave equation 
(|28|l . it becomes clear that orbits can also be seen as 4- 
dimcnsional guided waves by what could be described as 
a 4-dimensional optical fibre. Modes are to be expected 
in these waveguides and we shall say something about 
them later on. 



V. FERMAT'S PRINCIPLE IN 4 DIMENSIONS 

Fermat's principle applies to optics and states that 
the path followed by a light ray is the one that makes 
the travel time an extremum; usually it is the path that 
minimizes the time but in some cases a ray can follow a 
path of maximum or stationary time. These solutions are 
usually unstable, so one takes the view that light must 
follow the quickest path. In Eq. 1)31) |l we have defined a 
time interval associated with a 4-dimensional elementary 
displacement, which allows us to determine, by integra- 
tion, a travel time associated with displacements of any 
size along a given 4-dimensional path. We can then ex- 
tend Fermat's principle to 4D and impose an extremum 
requirement in order to select a privileged path between 
any two 4D points. Taking the square root to Eq. H30|l 



dt 



900 



dx'^dx^ . 



(41) 



Integrating between two points Pi and P2 

'■P2 I z rP2 



t 



Pi 



'^dx'dxi 
900 



Pi 



-^x'xJdt. (42) 
500 



In order to evaluate the previous integral one must know 
the particular path linking the points by defining func- 
tions x^{t), allowing the replacement dx* = iMt. At this 
stage it is useful to define a Lagrangian 



L = 



9i3 



25, 



00 



(43) 



(39) The time integral can then be written 



Since we have 54 normal to the other vectors we can ap- 
ply metric conversion and write the equivalent quadratic 
form for GTR 



dr' = [ — \ - ['-^ 



n4 

[dr^ +r^{de^ +81-0^ Od(p^)] 



(40) 



As we stated in the introduction, the usual form of 
Schwarzschild's metric is given by Eq. (P) but a more 
interesting, isotropic form is the one in Eq. Q. The 
latter can be compared to Eq. (|40|l allowing the deriva- 
tion of the refractive indices in Eqs. (0J). These refrac- 
tive indices provide a 4D0 Euclidean space equivalent 



t= / V2Ldt. 
'Pi 



(44) 



Time has to remain stationary against any small 
change of path; therefore we envisage a slightly distorted 
path defined by functions x^{t) + ex^{t), where e is arbi- 
trarily small and x*(^) ^.re functions that specify distor- 
tion. Since the distortion must not affect the end points, 
the distortion functions must vanish at those points. The 
time integral will now be a function of e and we require 
that 



dtje) 
de 



= 0. 



(45) 
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Now, the Lagrangian (|43|l is a function of , through gap 
and also an exphcit function of i*. AUowing for a path 
change, through e makes t in Eq. H44(l a function of e 



tie) 



P2 



^j2L{x^ + ex' + + ex*) dt. 



(46) 



Pi 



This can now be derived with respect to e 
di(e) 



de 



e=0 



1 



Pi 



2i 



dx^ dx^ ^ 



e=0 



Note that the first term on the rhs can be written 
1 dL 



P2 
Pi 



rx'dt 



/2Ldx 

This can be integrated by parts 



Pi dx^ 



(47) 



(48) 



Pi 



d{V2L) 
dx^ 




The first term on the second member is zero because 
vanishes for the end points; replacing in Eq. H47(l 



(50) 



dt{e) 






j-P-2 






de 


e=0 


x/2. 


IPl 


At \ 





1 dL 

L dx^ 



x'dt. 



The rhs must be zero for arbitrary distortion functions x*, 
so we conclude that the following set of four simultaneous 
equations must be verified 



d 
di 



1 dL 
VLdi" 



1 dL 
TLd^' 



(51) 



these are called the Euler-Lagrange equations. 

Consideration of Eqs. H31() and (|34|l allows us to con- 
clude that the Lagrangian defined by (|43|l can also be 
written as L = and must always equal 1/2. From 
the Lagrangian one defines immediately the conjugate 
momenta 



dL 

dx^ 



500 



(52) 



Notice the use of the lower index (w^) to represent mo- 
menta while velocity components have an upper index 
(w*). The conjugate momenta are the components of the 
conjugate momentum vector 



9 Vj 
V-900 



(53) 



and from Eq. © 



V-900V = 9'vi = 9'9iji^ = 9ji^- 



(54) 



The conjugate momentum and velocity are the same but 
their components are referred to the reciprocal and re- 
fractive index frames, respectively. [T^ Notice also that 
by virtue of Eq. H22|) it is also 



El 
Po 



(55) 



The Euler-Lagrange equations (|51|l can now be given 
a simpler form 



Vi 



d.,L. 



(56) 



This set of four equations defines trajectories of minimum 
time in 4D0 space as long as the frame vectors ga are 
known everywhere, independently of the fact that they 
may or may not be referred to the orthonormed frame 
via a refractive index. By definition these trajectories 
are the geodesies of 4D0 space, spanned by frame vectors 
5»/y-5oo, with metric tensor -5y /5oo- 

Following an exactly similar procedure we can find tra- 
jectories which extremize proper time, defined by taking 
the positive square root of Eq. (|34|l . The Lagrangian is 
now defined by 



c^-lii^x^'xr 
2 544 

Consequently the conjugate momenta are 
dxf^ g44 



(57) 



(58) 



From Eq. H22|) we have i'^ ~ p^/p4; the associated Euler- 
Lagrange equations are 



df^C. 



(59) 



"These are, by definition, spacelike geodesies of GTR 
with metric tensor —g^i^/ga and we have thus defined 
a method for one-to-one geodesic mapping between 4D0 
and spacelike GTR. Recalling the conditions for this 
mapping to be valid, all the frame vectors must be in- 
dependent of both t and r and go and 34 must be normal 
to the other 3 frame vectors. In tensor terms, all the gaf3 
must be independent from t and r and goi = g^4 = 0." 



VI. THE SOURCES OF REFRACTIVE INDEX 

The set of 4 equations (|56|l defines the geodesies of 4D0 
space; particularly in cases where there is a refractive 
index, it defines trajectories of minimum time but does 
not tell us anything about what produces the refractive 
index in the first place. Similarly the set of equations 
defines the geodesies of GTR space without telling 
us what shapes space. In order to analyse this question 



7 



we must return to the general case of a refractive frame 
without other impositions besides the existence of a 
refractive index. 

Considering the momentum vector 



(60) 



with Ug'^n ^j = 5^, we will now take its time derivative. 
Using Eq. (|D4l) 



p = x-(T)p) = x-G. 



(61) 



By a suitable choice of coordinates we can always have 
5" = (7°. We can then invoke the fact that for an elemen- 
tary particle in flat space the momentum vector com- 
ponents can be associated with the concepts of energy, 
3D momentum and rest mass as p — Ea^ + P + raa'^ 
(see d and Sec. IVITTl l If this consequence is ex- 
tended to curved space and to mass distributions, we 
write p = Ea^ + P + rng'^ , where now E is energy den- 
sity, p = Prnd^"" is 3D momentum density and m is mass 
density. The previous equation then becomes 



Ea'' 



mg 



x-G. 



(62) 



When the Laplacian is applied to the momentum vec- 
tor the result is still necessarily a vector 



D^p = S. 



(63) 



Vector S is called the sources vector and can be expanded 
into 25 terms as 



S = {B^n\)a[,p" = S^app"; 



(64) 



where = g"^pf^. Tensor S^p contains the coefficients 
of the sources vector and we call it the sources tensor. 
The sources tensor influences the shape of geodesies as 
we shall see in one particularly important situation. One 
important consequence that we don't pursue here is that 
by zeroing the sources vector one obtains the wave equa- 
tion D^p = 0, which accepts gravitational wave solutions. 

If a'^ is normal to the other frame vectors we can write 
p — E(a'^ + v) in the reciprocal frame, with v a unit 
vector OT p = E{—ao + v) in the direct frame. Equation 
H61|l can then be given the form 



E{a° + v) + Ev ^(70 + vG. 



(65) 



Since G can have scalar and bivector components, the 
scalar part must be responsible for the energy change, 
while the bivector part rotates the velocity v. The bivec- 
tor part of G is generated by DAp, which allows a sim- 
plification of the previous equation to 



V = v-(I)Av), 



(66) 



if the frame vectors are independent of t. This equation is 
exactly equivalent to the set of Euler-Lagrange equations 
H55|l but it was derived in a way which tells us when to 
expect geodesic movement or free fall. 



We will now investigate spherically symmetric solu- 
tions in isotropic conditions defined by Eq. this 
means that the refractive index can be expressed as func- 
tions of r. The vector derivative in spherical coordinates 
is of course 



D 



1 



(Jj-dr 



1 



credo 



1 



'Utdt ^ Grdr- 
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(67) 



The Laplacian is the inner product of D with itself but 
the frame vectors' derivatives must be considered; all the 
derivatives with respect to r are zero and the others are 



decre = -cTr, d^crg = cos^ct^, 

dga^p — 0, dipfftp — — sin 6 ar — cos 9 ag. 

After evaluation the curved Laplacian becomes 



(68) 



D^ 



■ dr 



dgg + 



+ - 



cot 9 „ csc^ 9 



■dg + 



a 

2 "^w 



(69) 



drr. 



The search for solutions of Eq. (|63|) must necessarily 
start with vanishing second member, a zero sources situ- 
ation, which one would implicitly assign to vacuum; this 
is a wrong assumption as we will show. Zeroing the sec- 
ond member implies that the Laplacian of both Ur and 
n4 must be zero; considering that they are functions of r 
we get the following equation for 



n H = 0, 



(70) 



with general solution — bexp(a/r). It is legitimate to 
make 6 = 1 because the refractive index must be unity 
at infinity. Using this solution in Eq. ()69|l the Laplacian 

becomes 



D^ = 



C-"/'- (^drr + ldr + ^dr + ^ Ogg + (71) 

cot 9 „ csc^ 9 „ \ „ 1 



CSC' t/ 



5, 



'tt 



(rn)' 



■ drr] 



which produces the solution ~ n^. So space must 
be truly isotropic and not relaxed isotropic as we had 
allowed. The solution we have found for the refractive 
index components in isotropic space can correctly model 
Newton dynamics, which led the author to adhere to it 
for some time 10] . However if inserted into Eq. 1)34(1 
this solution produces a GTR metric which is verifiably 
in disagreement with observations; consequently it has 
purely geometric significance. 

The inadequacy of the isotropic solution found above 
for relativistic predictions deserves some thought, so that 
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we can search for solutions guided by the results that are 
expected to have physical significance. In the physical 
world we are never in a situation of zero sources because 
the shape of space or the existence of a refractive index 
must always be tested with a test particle. A test particle 
is an abstraction corresponding to a point mass consid- 
ered so small as to have no influence on the shape of 
space; in reality a point particle is a black hole in GTR, 
although this fact is always overlooked. A test particle 
must be seen as source of refractive index itself and its 
influence on the shape of space should not be neglected 
in any circumstances. If this is the case the solutions for 
vanishing sources vector may have only geometric mean- 
ing, with no connection to physical reality. 

The question is then what should go into the second 
member of Eq. Ht)3|l in order to find physically meaning- 
ful solutions. If we are testing gravity we must assume 
some mass density to suffer gravitational influence; this is 
what is usually designated as non-interacting dust, mean- 
ing that some continuous distribution of non-interacting 
particles follows the geodesies of space. Mass density is 
expected to be associated with 5*^4; on the other hand 
we are assuming that this mass density is very small and 
so we use flat space Laplacian to evaluate it. We conse- 
quently make an ad hoc proposal for the sources vector 
in the second member of Eq. (|63|l 



S = — V^n4(T4. 



Equation H63() becomes 



-V n4(T4; 



(72) 



(73) 



as a result the equation for remains unchanged but 
the equation for becomes 



2n' 



2n' 



(74) 



When is given the exponential form found above, 
the solution is 714 = ^/n^. This can now be entered into 
Eq. (|34|l and the coefficients can be expanded in series 
and compared to Schwarzschild's for the determination 
of parameter a. The final solution, for a stationary mass 
M is 



„2M/r- 



714 



(75) 



The equivalent GTR space is characterized by the 
quadratic form 



(76) 



Expanding in series of M/r the coefficients of this metric 
one would find that the lower order terms are exactly 
the same as for Schwarzschild's and so the predictions of 
the metrics are indistinguishable for small values of the 
expansion variable. Montanus [12| | arrives at the same 
solutions with a different reasoning; the same metric is 
also due to Yilmaz [ll[ll|i3. 



Equation H73|l can be interpreted in physical terms as 
containing the essence of gravitation. When solved for 
spherically symmetric solutions, as we have done, the first 
member provides the definition of a stationary gravita- 
tional mass as the factor M appearing in the exponent 
and the second member defines inertial mass as V^n4. 
Gravitational mass is defined with recourse to some par- 
ticle which undergoes gravitational influence and is ani- 
mated with velocity v and inertial mass cannot be defined 
without some field acting upon it. Complete investi- 
gation of the sources tensor elements and their relation 
to physical quantities is not yet done; it is believed that 
16 terms of this tensor have strong links with homolo- 
gous elements of stress tensor in GTR, while the others 
are related to electromagnetic field. 



VII. ELECTROMAGNETISM IN 5D 
SPACETIME 

Maxwell's equations can easily be written in the form 
of Eq. H63|) if we don't impose the condition that 174 should 
remain normal the other frame vectors; as we have seen 
in section Hvl this has the consequence that there will be 
no GTR equivalent to the equations formulated in 4D0. 

We will consider the non-orthonormed reciprocal frame 
defined by 
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(77) 



where q and m are charge and mass densities, respec- 
tively, and A = A^a^^ is the electromagnetic vector po- 
tential, assumed to be a function of coordinates t and 
but independent of r. The associated direct frame has 
vectors 



ffAi = O-^ A^Gi, 

m 



54 = 0-4; 



(78) 



and one can easily verify that Eq. (0 is obeyed. The 
momentum vector in the reciprocal frame is p ~ Ea^ -\- 
PmO'™ + qAfj_a'^ + ma'^ and G in the second member of Eq. 
(|61|l is G = qDA. We will assume D A to be zero, as one 
usually does in electromagnetism; also D can be replaced 
by ^ because the vector potential does not depend on 
T. It is convenient to define the Faraday bivector F = 
f^A, similarly to what is done in Ref. ^]; the dynamics 
equation then becomes 



and rearranging 



p + qA = qx-F; 



p = qx-F — qA. 



(79) 



(80) 



The first term in the second member is the Lorentz force 
and the second term is due to the radiation of an accel- 
erated charge. 

Recalling the wave displacement vector Eq. I|D1(I we 
have now 

q 



dx = a^dx" 



m 



■ A^aidx^. 



(81) 
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This corresponds to a refractive index tensor whose non- 
zero terms are 



71 Q. — 1, n — Afj,. 



(82) 



According to Eq. I|64|l the sources tensor has aU terms 
null except for the following 



1 T^2 



m 



(83) 



where D is the covariant derivative given by 

D = g'^da = a^df, + (rx^ + ^ A^a'^)d4. (84) 

m 

We can then define the current vector J verifying 

AV^A = ^VF = J, (85) 

where 



(86) 



Please refer to 6, Chap. 7] or to ^8, Part 2] to see how 
these equations generate classical electromagnetism, par- 
ticularly how setting the current vector to zero generates 
electromagnetic waves. 



VIII. MONOGENIC FUNCTIONS AND 
QUANTUM MECHANICS 



Dirac equation has been derived from the 5- 
dimensional monogenic condition in previous works jj, 
ITTij : the motivation for returning to the subject here is 
the correction of the electro-dynamics equation, which is 
incorrect in the earlier paper and absent in the later 
one Q . Because we are working in geometric algebra, our 
quantum mechanics equations will inherit that character 
but the isomorphism between the geometric algebra of 
5D spacetime, 6*4.1, and complex algebra of 4 * 4 matri- 
ces, Af(4, C), ensures that they can be translated into 
the more usual Dirac matrix formalism. The equivalence 
between the two formulations has been amply demon- 
strated in the two references above. 

Recalling the monogenic condition (|19|l . we will now 
expand it into 3 terms 



(ct°(9o + CT™<9m + (T^di)ij = 0. 



(87) 



We have already established that this equation accepts 
solutions in the form of Eq. H21|l and we use that to eval- 
uate the derivative with respect to x'^ 



(88) 



If the equation is multiplied by a* on the left, the first 4 
terms on the first member acquire bivector factors of the 
form (T^^, the first of which, cr^" squares to unity, while 



the other 3 square to minus unity. Those bivectors be- 
long to the even sub- algebra of 6*4^1, which is isomorphic 
the the algebra of Minkowski spacetime, as we have al- 
ready stated. It is perfectly legitimate to replace the said 
bivector factors by Dirac matrices, as was demonstrated 
in the above cited references. We can then rewrite the 
monogenic condition as 



(89) 



which can be immediately recognized as Dirac's equation 
if P4 is assigned to the particle's rest mass. The mono- 
genic function given by Eq. l(T^ can then be given the 
usual physical interpretation of a Dirac spinor 



i(-Ei+p-x+mr) 



(90) 



where E is energy, p is 3-dimensional momentum and m 
is rest mass. 

In order to separate left and right spinor components 
we use a technique adapted from Ref. 0. We choose 
an arbitrary base element which squares to identity, for 
instance (T4, with which we form the two idempotents 
(1 -I- cr4)/2 and (1 — tT4)/2. The name idempotents 
means that they reproduce themselves when squared. 
These idempotents absorb any cr4 factor; as can be easily 
checked (1 + 0^)0^ = {1 + (T/^) and (1 — cr4)cr4 — — (1 — (74). 
Obviously we can decompose the wavefunction ip as 

^ = ^i±^ + ^i^ = ^,+ + ^_. (91) 

This apparently trivial decomposition produces some sur- 
prising results due to the following relations 



e'«(l + a4) = 



(cose* + isin6')(l + 0-4) 

(lcos6' + icr4sin6')(l +(74) (92) 



and similarly 



(93) 



We could have chosen other idempotents, which would 
produce similar results. The available idempotents gen- 
erate an SU (4) group and it has been argued that they 
may be related to different elementary particles. 'll*| 

Electrodynamics can now be implemented in the the 
same way used in Sec. IVIII to implement classical elec- 
tromagnetism. The monogenic condition must now be 
established with the covariant derivative given by Eq. 

(|H3 

(T^a^V + ('^'' + ^ ^M^^) diijj = 0. (94) 

Multiplying on the left by a'^ and taking d/^ilj = irml! 



[7^(5^ + igA^)+im]V = 0. 



(95) 



This equation can be compared to what is found in any 
quantum mechanics textbook.. 



It is now adequate to say a few words about quantiza- 
tion, which is inherent to 5D monogenic functions. We 
have already seen that these functions are 4-dimensional 
waves, that is, they have 3-dimensional wavefronts nor- 
mal to the direction of propagation. Whenever the re- 
fractive index distribution traps one of these waves a 4- 
dimensional waveguide is produced, which has its own 
allowed propagating modes. In the particular case of a 
central potential, be it an atom's or a galaxy's nucleus, 
we expect spherical harmonic modes, which produce the 
well known electron orbitals in the atom and have un- 
known manifestations in a galaxy. 
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FIG. 1: Indices in the range {0, 4} will be denoted with 
Greek letters a, /3, 7. Indices in the range {0, 3} will also re- 
ceive Greek letters but chosen from fi, v, ^. For indices in the 
range {1, 4} we will use Latin letters k and finally for in- 
dices in the range {1,3} we will use also Latin letters chosen 
from m, n, o. 



IX. CONCLUSION 

Every physicist dreams of finding a unified formula- 
tion for the fundamental laws of physics. It is usually 
accepted that in order to achieve such objective a new 
paradigm is needed, meaning that one must surely step 
back from accepted physics principles and start afresh 
from new simpler ones. Ideally one should have a small 
set of principles, valid for all areas of physics, and all the 
important relations should flow naturally from mathe- 
matical reasoning. 

In this paper we extend a proposal previously made 
in that direction, that one should accept 5-dimensional 
spacetime as the adequate space to formulate the laws 
of physics and introduce in this space the condition of 
monogeneity. We had shown in another work that this 
condition is sufficient to arrive simultaneously at special 
relativity and the free particle Dirac equation; here we 
show that by generalizing the monogenic condition to 
bent spaces one is able to obtain relativistic dynamics not 
entirely coincident with GTR but also electrodynamics, 
both classical and quantized. 

Maxwell's equations were also derived from the mono- 
genic condition and were unified to the equations respon- 
sible for gravitational dynamics. The procedure is not yet 
entirely satisfactory, in the sense that an ad hoc proposal 
had to be made in respect to inertial mass; in future work 
we hope to find a suitable formulation for the derivation 
of curvature from first principles. 



APPENDIX B: NON-DIMENSIONAL UNITS 

The interpretation of t and r as time coordinates im- 
plies the use of a scale parameter which is naturally cho- 
sen as the vacuum speed of light c. We don't need to 
include this constant in our equations because we can 
always recover time intervals, if needed, introducing the 
speed of light at a later stage. We can even go a step fur- 
ther and eliminate all units from our equations so that 
they become pure number equations; in this way we will 
avoid cumbersome constants whenever coordinates have 
to appear as arguments of exponentials or trigonomet- 
ric functions. We note that, at least for the macroscopic 
world, physical units can all be reduced to four funda- 
mental ones; we can, for instance, choose length, time, 
mass and electric charge as fundamental, as we could 
just as well have chosen others. Measurements are then 
made by comparison with standards; of course we need 
four standards, one for each fundamental unit. But now 
note that there are four fundamental constants: Planck 
constant {Ti), gravitational constant (G), speed of light 
in vacuum (c) and proton electric charge (e), with which 
we can build four standards for the fundamental units. 
Table lists the standards of this units' system, fre- 

TABLE I: Standards for non-dimensional units' system 



Length 



Time 



Mass 



Charge 



APPENDIX A: INDEXING CONVENTIONS 

In this section we establish the indexing conventions 
used in the paper. We deal with 5-dimensional space 
but we are also interested in two of its 4-dimensional 
subspaces and one 3-dimensional subspace; ideally our 
choice of indices should clearly identify their ranges in 
order to avoid the need to specify the latter in every 
equation. The diagram in Fig.^shows the index naming 
convention used in this paper; Einstein's summation con- 
vention will be adopted as well as the compact notation 
for partial derivatives da = d/dx"^. 



quently called Planck units, which the authors prefer to 
designate by non-dimensional units. In this system all 
the fundamental constants, Ti, G, c, e, become unity, a 
particle's Compton frequency, defined hy v = mc^ /h, be- 
comes equal to the particle's mass and the frequent term 
GM/{(?r) is simplified to M/r. We can, in fact, take 
all measures to be non-dimensional, since the standards 
are defined with recourse to universal constants; this will 
be our posture. Geometry and physics become relations 
between pure numbers, vectors, bivectors, etc. and the 
geometric concept of distance is needed only for graphical 
representation. 
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APPENDIX C: SOME COMPLEMENTS OF 
GEOMETRIC ALGEBRA 

In this section we expand the concepts given in Sec.lTTl 
introducing some useful relations and definitions. Start- 
ing with the basis elements that square to unity Eq. ((HJ, 
repeated here, 



(CI) 



it is easy to verify any of the above equations; sup- 
pose we want to check that (croy)^ = 1. Start by ex- 
panding the square and remove the compact notation 
((Toy)^ — '^o'^i'^j'^oo'i'^j , then swap the last aj twice to 
bring it next to its homonymous; each swap changes the 
sign, so an even number of swaps preserves the sign: 
(foy)^ — o'aO'iio'j)'^aQai. From the third equation 
we know that the squared vector is unity and we get suc- 



cessively (cTOij) = (ToCTiaoCTi = -{aoyia^y 



using the first equation ^ we get finally (croy)^ = 1 as 
desired. 

The remaining basis elements square to —1 as can be 
verified in a similar manner, Eq. (jT)): 

{aof - = {iaif = i^ = -1. (C2) 

Note that the pseudoscalar i commutes with all the other 
basis elements while being a square root of —1; this makes 
it a very special element which can play the role of the 
scalar imaginary in complex algebra. 

We can now address the geometric product of any two 
vectors a = a'^aa and b — b^a/s making use of the dis- 
tributive property 



ab 



-a"b 



a'b' 



a/3; 



(C3) 



and we notice it can be decomposed into a symmetric 
part, a scalar called the inner or interior product, and 
an anti-symmetric part, a bivector called the outer or 
exterior product. 



ab^a-b + aAb, ba = a-b — aAb. 



(C4) 



Reversing the definition one can write inner and outer 
products as 



1 



1 



a-b = — (ab + ba), aAb — — (ab — ba). 



(C5) 



The inner product is the same as the usual "dot prod- 
uct," the only difference being in the negative sign of 
the aobo term; this is to be expected and is similar to 
what one finds in special relativity. The outer prod- 
uct represents an oriented area; in Euclidean 3-space 
it can be linked to the "cross product" by the relation 
cross(a, b) — — (7i23aAb; here we introduced bold char- 
acters for 3-dimensional vectors and avoided defining a 
symbol for the cross product because we will not use it 



again. We also used the convention that interior and ex- 
terior products take precedence over geometric product 
in an expression. 

When a vector is operated with a multivector the inner 
product reduces the grade of each element by one unit 
and the outer product increases the grade by one. We 
will generalize the definition of inner and outer products 
below; under this generalized definition the inner product 
between a vector and a scalar produces a vector. Given 
a multivector a we refer to its grade-r part by writing 
<a>r; the scalar or grade zero part is simply designated 
as < a >. By operating a vector with itself we obtain a 
scalar equal to the square of the vector's length 



a = aa — a-a -\r aAa ^ a-a. 



(C6) 



The definitions of inner and outer products can be ex- 
tended to general multivectors 



a.l3 



aAb = ^(<a>„ <6>/3)„+^. 

a,l3 



(C7) 
(C8) 



Two other useful products are the scalar product, de- 
noted as <ab> and commutator product, defined by 



X b = ab — ba. 



(C9) 



In mixed product expressions we will use the convention 
that inner and outer products take precedence over geo- 
metric products. 

We will encounter exponentials with multivector expo- 
nents; two particular cases of exponentiation are specially 
important. If u is such that = — 1 and is a scalar 

= — + — -...{= cos9} 
2! 4! ^ ' 

-^^6*-^— + ...{= u sin 6*} (CIO) 
— cos 6* -I- u sin 6*. 

Conversely if h is such that — 1 



q3 m 



e- = ^ + hO+^ + h- + - + ... 

+he + h— + ...{=hsin\ie] (cii) 

= cosh 6 + h sinh 0. 

The exponential of bivectors is useful for defining rota- 
tions; a rotation of vector a by angle 9 on the (T12 plane 
is performed by 



RaR; 



(C12) 
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the tilde denotes reversion and reverses the order of all 
products. As a check we make a — ai 



cos - 
2 



0-12 sm - 



* I cos - + (Ti2 sm - ) (C13) 

= cos 9(7i + sin 6a2 . 

Similarly, if we had made a — a2, the result would have 
been — sin 9ai + cos 9a2 ■ 

If we use B to represent a bivector whose plane is nor- 
mal to ao and define its norm by \B\ — {BBY^^, a general 
rotation in 4-space is represented by the rotor 



= = cos 



R = e 



\B\ 



B 



The rotation angle is and the rotation plane is defined 
by i?. A rotor is defined as a unitary even multivector 
(a multivector with even grade components only) which 
squares to unity; we are particularly interested in rotors 
with bivector components. It is more general to define 
a rotation by a plane (bivector) then by an axis (vector) 
because the latter only works in 3D while the former is 
applicable in any dimension. When the plane of bivector 
B contains ctq, a similar operation does not produce a 
rotation but produces a boost instead. Take for instance 
B = (7010/2 and define the transformation operator T = 
exp(i?); a transformation of the basis vector ctq produces 

a' = f aoT = e-"«i''/Voe"°i''/2 



cosh cToi sinh - ) ao 



* ^cosh — + (jgi sinh — 
— cosh 9ao + sinh 9ai . 



(C15) 



APPENDIX D: TIME DERIVATIVE OF A 
4-DIMENSIONAL VECTOR 

If there is a refractive index the wave displacement 
vector can be written as 



da; — gadx" — n^acrfsdx' 



(Dl) 



Because this vector is nilpotent, by virtue of Eq. (|29|l . the 
five coordinates are not independent and we can divide 
both members by da;" = dt defining the nilpotent vector 



i — .90 + 5i2;* = n"oCa + n^iapx^ 



(D2) 



Suppose we have a 5D vector a — aaa'^ and we want 
to find its time derivative along a path parameterized by 
t, that is all the a;' are functions of t. We can write 

a = dpa'^x^aa, (D3) 
where naturally x'^ — 1. Remembering the definition of 
covariant derivative 1)16(1 and Eq. (|D2p we can modify this 
equation to 



x^gp-g'^dpa'^aa = x-{Da). 



(D4) 



We have expressed vector a in terms of the orthonormed 
frame in order to avoid vector derivatives but the result 
must be independent of the chosen frame. 

This procedure has an obvious dual, which we arrive 
at by defining 

x = 5^x^+54. (D5) 
The proper time derivative of vector a is then 

a = x-(Da). (D6) 
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